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Neuroimaging studies of post-traumatic stress disorder (PTSD)-related memory impairments have consis- 
tently implicated abnormal activities in the frontal and parietal lobes. However, most studies have used 
block designs and could not dissociate the multiple phases of working memory. In this study, the involve- 
ment of the prefrontal cortex in working memory phases was assessed among veterans with PTSD and age-/ 
gender-matched healthy controls. Multichannel functional near infrared spectroscopy (fNIRS) was utilized 
to measure prefrontal cortex hemodynamic activations during memory of neutral (i.e., not trauma-related) 
forward and backward digit span tasks. An event-related experimental design was utilized to dissociate the 
different phases (i.e., encoding, maintenance and retrieval) of working memory. The healthy controls showed 
robust hemodynamic activations during the encoding and retrieval processes. In contrast, the veterans with 
PTSD were found to have activations during the encoding process, but followed by distinct deactivations 
during the retrieval process. The PTSD participants, but not the controls, appeared to suppress prefrontal 
activity during memory retrieval. This deactivation was more pronounced in the right dorsolateral prefrontal 
cortex during the retrieval phase. These deactivations in PTSD patients might implicate an active inhibition 
of dorsolateral prefrontal neural activity during retrieval of working memory. 

© 2014 The Authors. Published by Elsevier Inc. 
This is an open access article under the CC BY-NC-ND license 
(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Post-traumatic stress disorder (PTSD) is a neuropsychiatric dis- 
order that can develop after exposure to any traumatic event. The 
estimated lifetime prevalence of PTSD is 7.8% among adult Americans 
(Kessler et al., 1995). Military combat is the most common trauma 
among men with PTSD. Combat-related PTSD is found in 9%-25% of 
war-zone veterans (Kessler et al, 1995; Dohrenwend et al, 2006; 
Erickson et al, 2001 ; Hoge et al, 2004) and is often persistent and co- 
morbid with other mental disorders (Arbanas, 2010), even after the 
veterans return to their civilian lives. 

Numerous neuropsychological studies have identified cognitive 
dysfunctions associated with PTSD, such as memory impairments, 
attention deficits, and learning disabilities (Bremner et al, 1993; 
Vasterling et al, 1998, 2002; Brandes et al, 2002; Koenen et al, 
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2001 ; Roca et al, 2006). In recent years, various neuroimaging tech- 
niques have been used to investigate neural mechanisms underlying 
stress-related and trauma-related memory impairments. Such imag- 
ing techniques include event-related potentials (ERPs) (Galletly et al, 
2001; Weber et al, 2005; Veltmeyer et al, 2009), positron emission 
tomography (PET) (Shaw et al, 2002; Clark et al, 2003), functional 
magnetic resonance imaging (fMRI) (Elzinga et al, 2007; Moores et 
al, 2008; Shaw et al, 2009), and functional near infrared spectroscopy 
(fNIRS) (Matsuo et al, 2003). The most replicable findings from the 
neuroimaging studies suggest abnormal activities in the frontal and 
parietal cortices, two of the key cerebral regions consistently impli- 
cated in the neurobiology of PTSD (Bremner, 2002; Francati et al, 
2007; Aupperle et al, 2012). 

Working memory involves temporary brain storage and manip- 
ulation of information necessary for executing complex cognitive 
tasks. It is composed of multiple component processes or phases, 
such as encoding, maintenance, and retrieval of the given informa- 
tion (Baddeley, 1992, 1996). Neural correlates of these processes can 
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be temporally dissociated and identified by using event-related ex- 
perimental designs (Aguirre and D'Esposito, 1999). However, most 
studies of PTSD-related working memory using PET, fMRI, and fNIRS 
techniques have used block designs. Block designs do not allow reli- 
able dissociation of overlapping brain responses during multi-phase 
processes of working memory. Subtracting brain responses during 
blocked trials with and without specific memory processes are often 
challenged by the untenable assumption of pure insertion (Zarahn et 
al., 1999; D'Esposito et al, 1999). While ERP provides a direct mea- 
sure of event-related electrophysiological response of the brain, the 
source of a single evoked potential cannot be easily localized and iso- 
lated. To date, little is known on how temporal and spatial responses 
of the brain are altered by PTSD during the multi-phase processes of 
working memory. 

Functional near infrared spectroscopy (fNIRS) is a noninvasive 
neuroimaging technology that measures the cerebral hemodynamics 
and oxygenation with near infrared light (670-900 nm) (Ferrari et al., 
2012; Boas et al, 2004). The near infrared light can penetrate through 
the scalp and skull to reach a depth of several centimeters within 
the cortex (Dehghani et al., 2009; Cui et al, 201 1 ; Tian and Liu, 2014), 
while being partially absorbed by the oxygenated hemoglobin (Hb0 2 ) 
and deoxygenated hemoglobin (Hb) in the cerebral blood flow. By 
measuring changes of near infrared light intensity from the baseline 
at two or more wavelengths, the changes of HDO2 and Hb concentra- 
tions in the cerebral vasculature can be quantified (Cope et al, 1988). 
Compared with other neuroimaging modalities, this technology has 
an excellent temporal resolution (in ms) and also has such advantages 
as being portable, relatively inexpensive, and feasible to be used in 
natural environments without too much body confinement (Hoshi et 
al., 2007). Recent studies have demonstrated that fNIRS could be uti- 
lized effectively and successfully to study prefrontal cortex activation 
associated with the working memory processes in human participants 
(Moro et al, 2013; Sato et al., 2013). 

In this study, we used fNIRS to image the right and left prefrontal 
activations during working memory phases (encoding, maintenance, 
and retrieval) in veterans with PTSD and compared their prefrontal re- 
sponses with those in age- and gender- matched healthy controls. We 
chose the forward and backward digit span tasks because they have 
been intensively used to assess working memory in neuropsychologi- 
cal research (Gerton et al, 2004). Both tasks require immediate, serial 
recall of a list of digits that were presented to the participants. In a 
digit forward task, participants repeat the digits in the same order 
as they are presented; in a digit backward task, participants repeat 
the digits in reverse order. A digit backward task is believed to share 
many of the cognitive components of a digit forward task except that 
the digit backward task requires an additional manipulation compo- 
nent. We used an event-related experimental design to image and 
dissociate the multi-phase processes of working memory during the 
digit span tasks. In this way, we were able to reveal both temporal 
and spatial patterns of prefrontal responses that could be specifi- 
cally altered in patients with PTSD. The overall goal of this study was 
to demonstrate that fNIRS could be a portable and complementary 
neuroimaging tool to identify abnormal prefrontal activities in PTSD 
patients during working memory processes. 

2. Material and methods 

2.1. Participants 

2.1.1. Veterans with PTSD 

Combat-exposed male veterans who were preparing for or en- 
rolled in higher education were recruited via a supported education 
intervention clinical trial conducted by the mental health transla- 
tional research center of the School of Social Work at the University of 
Texas at Arlington. The veterans were referred for the fNIRS scan (to be 
described in the following Subsections 2.2 and 2.3) if they had medical 



records confirming a prior diagnosis with PTSD, and were experienc- 
ing clinically significant distress and functional impairment affecting 
their cognitive and related academic performance at the time of study. 
Confirmation of their PTSD diagnosis and comorbid conditions was 
conducted prior to the fNIRS scan by an experienced licensed clinician 
using the Structured Clinical Interview for DSM Disorders — research 
version (SCID-RV, http://www.scid4.org/info/refscid.html) (First et 
al, 2002) and scores on three self-report questionnaires: the PTSD 
checklist - military version (PCL-M) (Weathers et al, 1993; Wilkins 
et al, 2011), Mississippi scale for combat-related PTSD (Hyer et al, 
1991; Keane et al, 1988), and the PTSD subscale of the Minnesota 
Multiphasic Personality Inventory (Keane et al, 1984). At least one of 
the three self-report scores was required to be above its respective 
clinical cutoff value (i.e., PCL-M > 50, Mississippi scale > 1 30, or PTSD 
subscale > 30) for the veterans to be enrolled, besides their prior di- 
agnosis and current diagnostic evaluation confirmation of PTSD and 
comorbid conditions. 

A total of 17 veterans who met the eligibility criteria were suc- 
cessfully recruited and scanned by fNIRS. One veteran's results were 
dropped from the study database during data analysis because all of 
his trials of correct responses had motion artifacts. For the remaining 
16 veterans, current comorbid conditions included attention deficit 
hyperactivity disorder (ADHD) (n = 6), major depressive disorder 
(n = 5), alcohol dependence (n = 4), musculoskeletal pain (n = 3), 
insomnia (n = 3), history of blast exposure (n = 3), anxiety disor- 
der (n = 2), mild traumatic brain injury (mTBI) (n = 2), and learning 
disorder (n = 2). 

2. 1 .2. Healthy controls 

A total of 17 healthy male adults were also recruited from the 
local community of the University of Texas at Arlington and matched 
with the PTSD group on gender and age. Potential participants were 
screened initially with a questionnaire to exclude those who had 
history of any psychiatric illness, history of severe TBI, or intake of 
any medication at the time of study. Corresponding to the veteran 
who was excluded during data analysis, we also excluded the data 
from his healthy counterpart, yielding data from 16 healthy controls. 

All of the 32 remaining participants in data analysis were 
right-handed, with age ranging from 21 to 56 years, mean ± SD 
age = 29.4 ± 9.6 years. The study protocols for both veterans with 
PTSD and matched healthy controls were approved by the University 
of Texas at Arlington Institutional Review Board. Written informed 
consent was obtained from every participant prior to the fNIRS scan. 

2.2. Tasks and paradigm 

All participants sat comfortably before a computer, and completed 
a session of digit forward task (8 trials) and a session of digit backward 
task (8 trials) sequentially while they were scanned by fNIRS. As mod- 
ified from a previous study (Sun et al, 2005), each trial started from a 
string of six digits (randomly generated by the computer) that were 
presented in sequence on an LCD monitor, one digit per second. Then 
a blank screen was displayed for 10 s, during which the participants 
were instructed to look at the screen and covertly rehearse the given 
digits continuously. After this retention interval, a number pad was 
displayed on the monitor. The participants were instructed to recall 
the digits either in the same order (forward task) or reverse order 
(backward task) as they were presented, by clicking the on-screen 
number buttons. All participants were able to finish the forward or 
backward recall within 10 s. At last, an inter-trial interval of 10 s was 
given before the next trial. According to the study by Sun et al. (2005), 
a 10-s inter-trial interval would be sufficient for the hemodynamic 
signal to recover to the baseline level, while maintaining participants' 
patience on the task. 

Before each formal measurement session was started, all partici- 
pants were trained to practice a few trials. An experimenter observed 
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Fig. 1. Configuration and cortical position of the fNIRS probe: (A) configuration of the 
fNIRS probe. Red circles represent light sources, blue squares represent detectors, and 
green lines represent the nearest source-detector pairs (channels) to measure the brain 
activities. (B) Co-registered positions of the optodes on a standard brain atlas. The probe 
partially covered the anterior, dorsolateral, and ventrolateral prefrontal regions on the 
right and left hemispheres. The anatomical position of each channel on the brain atlas 
is reported in detail in Table 1. 

the course of practice to confirm that the participants understood the 
protocol and were able to complete the task correctly. 

It is known that for both digit forward and backward tasks, the 
task difficulty depended on the length of digit string (Gerton et al, 
2004). We selected six-digit tasks because they are at an intermediate 
level and can be completed by all of the control and PTSD participants 
without extreme difficulty. In addition, the participants were always 
instructed to complete the digit forward task first. Because the digit 
forward task was easier than the digit backward task, this sequence 
allowed the participants to adapt to the following more difficult task 
ssession. The accuracy of each participant's performance in each type 
of task was measured by the percentage of trials of correct retrieval 
relative to the total number of trials (=8) in each session. 

2.3. Functional near infrared spectroscopy 

2.3.1. Instrument 

A high-performance fNIRS system (Cephalogics LLC, Boston, MA) 
(Zeff et al, 2007) was used to acquire each participant's prefrontal 
hemodynamic activities during performance of the tasks. The sys- 
tem used light-emitting diodes (LEDs) at two wavelengths (750 and 
850 nm) as light sources and avalanche photodiodes (APDs) as de- 
tectors. The data sampling rate was 10.8 Hz. The configuration of 
the fNIRS probe is shown in Fig. 1A. It was composed of 12 light 
sources and 16 detectors (6 light sources and 8 detectors on each 
hemisphere). The probe was placed bilaterally and symmetrically on 
the participant's forehead. The bottom line of 6 light sources in the 
probe was just above the eyebrows, and its midpoint was 3.5 cm in 
distance from the nasion (see Fig. 1A). 

The fNIRS probe provided a total of 36 measurements (channels) 
when only the nearest source-detector pairs were considered (the 
nearest source-detector distance was 2.8 cm). Other larger source- 
detector distances were not included because their signals were too 
weak to be scientifically meaningful. The probe assembly was con- 
structed with low-weight optical fibers (TechEn Inc., Boston, MA) and 
thin polyethylene film to ensure participants' comfort during the ex- 
periment. 

2.3.2. Spatial registration 

To estimate the cortical regions that the fNIRS probe covered, 
a spatial registration procedure (Singh et al, 2005) was performed 
among six randomly selected participants: Once the fNIRS probe was 
set in place, the positions of optodes (i.e., light sources and detectors) 
along with five cranial landmarks (the nasion, inion, left and right 
preauricular points, and vertex) were measured using a PATRIOT mo- 
tion tracking system (Polhemus, Colchester, Vermont, USA). The cra- 
nial landmarks served as mediators to convert the real-world stereo- 
taxic coordinates of the optodes to the Montreal Neurological Institute 
(MNI) coordinates used in a standard brain MRI atlas (Tian et al, 2012) 



Table 1 



Registered positions of fNIRS measurement channels on the standard brain MRI atlas. 
The channel numbers are defined in Fig. 1A. 


Hemisphere 


Brodmann area 


Channel 


Right 


BA9 


13, 14, 15, 16, 17, 18 




BA10 


4, 5,6 




BA 44 


9,10 




BA45 


1,7,8 




BA46 


2,3, 11, 12 


Left 


BA 9 


31,32,33,34,35,36 




BA10 


22, 23, 24 




BA 44 


27,28 




BA45 


19,25,26 




BA46 


20,21,29,30 



based on the affine transformation. Fig. IB shows the registered op- 
tode positions (averaged over six participants) on the standard brain 
atlas. The probe partially covered the anterior, dorsolateral and ven- 
trolateral prefrontal regions, including Brodmann areas (BAs) 9, 10, 
44, 45, and 46. A detailed report is given in Table 1. Because the fNIRS 
probe was carefully placed on each individual participant's forehead 
by referring to the nasion and eyebrows, the registered optode posi- 
tions across individual participants were relatively consistent (posi- 
tional variations were about 5 mm) in comparison with the separation 
of two neighboring measurement channels (about 12.5 mm). 

2.3.3. Data screening and processing 

The fNIRS data from each measurement session was screened and 
processed using a publically available toolbox (Huppert et al, 2009). 
First, the raw data were visually inspected to exclude trials asso- 
ciated with significant data discontinuities (with a signal swing of 
15% or larger from the baseline intensities), which were usually at- 
tributed to motion artifacts during the experiment. Then the qualified 
data in optical density were low-pass filtered at a cut-off frequency of 
0.2 Hz to remove electronic noise and fast-oscillating cardiac waves, 
and high-pass filtered at a cut-off frequency of 0.01 Hz to remove 
possible slow baseline drift. Then, changes of hemoglobin concentra- 
tions relative to the baseline, A[Hb02] and A[Hb], were computed 
following the modified Beer-Lambert Law (Cope et al, 1988). At this 
step, we estimated the differential pathlength factor (DPF) to be 6.2 
at 750 nm and 5.8 at 850 nm, according to published data for adult 
heads (Essenpreis et al, 1993). Event-related temporal profiles were 
averaged over the remaining trials in each session to obtain individual 
hemodynamic responses (in A[Hb02] and A[Hb]) that were evoked 
by the digit forward and backward tasks. At this step, only the tri- 
als of correct retrievals, in which the participants' competence of the 
memory processes could be confirmed, were averaged. 

2.3.4. General linear model (GLM) analysis 

GLM analysis has been increasingly utilized to analyze fNIRS data 
over the last decade to identify cortical areas significantly stimulated 
by given tasks (Plichta et al, 2006; Ye et al, 2009; Tian and Liu, 2014; 
Hassanpour et al, 2014). In this work, we also employed GLM to ana- 
lyze channel-wise, event-related hemodynamic responses in the pre- 
frontal cortex stimulated by both digit forward and backward tasks. In 
particular, we used multiple regression analysis to account for multi- 
phase working memory performance. Similar to Sun et al. (2005), 
three memory phases were considered in each of the two tasks: stim- 
ulus (or encoding), delay (or active maintenance), and recall (or re- 
trieval). Accordingly, three regressors (encoding, maintenance, and 
retrieval) were generated after convolving three corresponding box- 
car functions with a hemodynamic response function (HRF) (Tian and 
Liu, 2014). Ideally, the HRF would be derived from the fNIRS sig- 
nals via an event-related experimental paradigm. To the best of our 
knowledge, however, such an HRF is not available for memory-evoked 
responses in the prefrontal cortex. Therefore, we used a standard HRF 
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Fig. 2. Schematic diagram of GLM analysis with three regressors that were generated 
to fit the hemodynamic responses evoked during the encoding, maintenance, and 
retrieval phases of working memory. 

derived from BOLD flvIRI as a surrogate (Glover, 1999). By fitting the 
three regressors to the channel-wise, event-related hemodynamic 
responses, the amplitudes (expressed by ^-values) of prefrontal ac- 
tivations or deactivations in responses to each memory phase were 
estimated by GLM (Tian and Liu, 2014). A schematic diagram of GLM 
analysis with three regressors is shown in Fig. 2. 

2.4. Statistical analysis 

For behavioral measures, the Wilcoxon rank-sum test was used 
to compare performance accuracies between the control and PTSD 
groups as well as between the forward and backward tasks. The 
Wilcoxon rank-sum test is a nonparametric alternative to the t test 
and is often used when the data cannot be assumed to be normally 
distributed. 

For hemodynamic measures, the channel-wise ^-values were ex- 
amined with the interaction model of analysis of variance (ANOVA). 
Three factors were considered: group (control and PTSD), task (for- 
ward and backward), and phase (encoding, maintenance, and re- 
trieval). If the main effect and/or interaction were found to be sig- 
nificant, separate t-tests were performed to identify the regions of 
interest (ROIs) that showed the significant difference. Both t-statistic 
values (expressed by t-values) and p-values were derived from the 
t-tests for each channel. 

Lastly, correlations between the behavioral and hemodynamic 
measures were tested based on linear regression analyses. These 
analyses determined if the hemodynamic measures could be ob- 
jective predictors of the participants' performance at the individual 
level. Here the hemodynamic changes were computed as the mean 
/3 -values within the respective ROIs identified above. 

2.5. Topography 

Topographic images of prefrontal activation were generated using 
EasyTopo, an optical topography toolbox we developed recently (Tian 
et al, 2012). EasyTopo is based on a standard brain MRI atlas and 
implements 2D angular interpolation of the channel-wise data in a 
spherical coordinate system. In this study, the channel-wise /3 -values 
derived from GLM analysis and t-values from subsequent statistical 
comparisons were interpolated to generate activation maps (/3-maps 
and t-maps) in each memory phase. 

3. Results 

3. 1 . Behavioral measures 

Table 2 reports the performance accuracies of control and PTSD 
groups for the digit forward and backward tasks, respectively. The 



healthy controls had higher accuracies than the veterans with PTSD in 
both tasks. The group difference was highly significant in the forward 
task (p < 0.01) and significant in the backward task (p = 0.04 < 
0.05). In terms of task differences, the healthy controls had reasonably 
higher accuracies in the forward task than those in the backward task 
(p < 0.01 ). In contrast, the veterans with PTSD had approximately 
equal accuracies in both tasks (p = 0.36 > 0.05). 

3.2. Hemodynamic responses 

For hemodynamic responses, ANOVA analysis revealed that the 
group effect (control vs. PTSD) was statistically significant in 28 of the 
total 36 channels (p < 0.01). The interaction between group (control 
vs. PTSD) and task (forward vs. backward) was also significant in 12 
of the total 36 channels (p < 0.05). Therefore, the group difference 
within each task and task difference within each group were identified 
in separate t-tests. 

3.2.1. Hb02 changes induced by forward and backward tasks in control 
group 

For the healthy controls, the digit forward task evoked robust Hb0 2 
changes, while the Hb changes were relatively low and insignificant 
for all three phases, as shown in Fig. 3A. Specifically, the task-evoked 
Hb0 2 time course clearly showed two peaks corresponding to the 
encoding and retrieval phases (Fig. 3 A), which are in accordance with 
a previous fMRI report using a similar stimulation protocol (Sun et al, 
2005). Furthermore, the respective topographic images (first row in 
Fig. 4) showed that the regions of activation were primarily located 
in the dorsolateral prefrontal cortex (DLPFC, BAs 9 and 46), and the 
left hemisphere had greater activation than the right hemisphere. 

In the control group, the Hb02 changes evoked by the digit back- 
ward task (Fig. 3B) had a much higher peak (D.58 \xM or microM) 
in the retrieval phase than those (OA \iM or microM) evoked by the 
digit forward task (paired t-test, p < 0.05). The digit backward task 
also evoked a subtle peak of Hb0 2 changes in the maintenance phase. 
According to the topographic images (second row in Fig. 4), the main 
regions of activation evoked by the digit backward task were also lo- 
cated in the DLPFC. The predominant difference between the forward 
and backward tasks was seen in both the DLPFC (BAs 9 and 46) and the 
anterior prefrontal cortex (APFC, BA 10) in the retrieval phase (third 
and fourth rows in Fig. 4). 

3.2.2. Hb0 2 changes induced by forward and backward tasks in the PTSD 
group 

For the veterans with PTSD, all Hb changes evoked by the forward 
and backward tasks were low and insignificant (Fig. 5A and B). The 
Hb0 2 changes evoked by the two types of tasks were very similar to 
each other: a subtle shoulder in the encoding phase before reaching 
a predominant peak in the maintenance phase, followed by a dis- 
tinct dip in the retrieval phase. In the forward task, the dip of Hb0 2 
changes in the retrieval phase reached below the baseline. The to- 
pographic images (first and second rows in Fig. 6) of the PTSD group 
showed a focused response region in the left DLPFC during the encod- 
ing and maintenance phases, and a distinct deactivation towards the 
APFC (BA 10) in the retrieval phase especially when completing the 
forward task. Additionally, the Hb0 2 changes evoked by the forward 
and backward tasks were similar to each other, with only a narrow 
region in APFC (BA 10) showing a significant task difference in the 
retrieval phase (third and fourth rows in Fig. 6). 

3.2.3. Control group vs. PTSD group 

For the digit forward task, since both the control and PTSD groups 
had activations during encoding and maintenance, between-group 
comparisons did not identify any channel with significance differ- 
ences (p < 0.05) in these two phases. In the retrieval phase, because 
the control group had a distinct activation and the PTSD group showed 
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Table 2 

Task-dependent performance accuracy of control and PTSD groups. 





Control group 




PTSD group 






Forward 


Backward 


Forward 


Backward 


Accuracy (%) 


95.0 ± 7.9 


78.3 ± 20.0 


70.3 ± 25.4 


62.5 ± 27.8 




Fig. 3. Mean task-evoked prefrontal hemodynamic responses (channel 30 in Fig. 1 A) in the control group: (A) mean hemodynamic responses to the digit forward task, in which the 
mean Hb0 2 changes show two robust peaks in the encoding and retrieval phases; (B) mean hemodynamic responses to the digit backward task, in which the mean Hb0 2 changes 
show higher magnitudes in the maintenance and retrieval phases compared with those in the digit forward task. In both panels, the solid lines represent the mean time courses 
of Hb0 2 , the dotted lines represent the mean time courses of Hb, the shaded regions represent the standard errors, the * symbols indicate the period of significant Hb0 2 changes 
(one-sample t-test, p < 0.01 ) from the baseline, and the gray bars in the bottom indicate the three phases of the task. It is noted that the actual retrieval time was slightly variable 
among participants. Therefore the retrieval phases labeled in both panels are schematic and approximate. 



a distinct deactivation, a large portion of the channels showed highly 
significant differences between the two participant groups (p < 0.01 ). 
According to the topographic images in Fig. 7, two ROIs with signif- 
icant group differences were identified in the retrieval phase: 1) a 
broader region in the right prefrontal cortex, which involved mainly 
the right middle frontal gyrus, and 2) a narrower region in the left 
prefrontal cortex, which was mainly located within the left superior 
frontal gyrus. 

For the digit backward task, there were also a large portion of 
the channels having highly significant between-group differences in 
the retrieval phase (p < 0.01 ). As shown in Fig. 8, two identified ROIs 
included a broader region in the right prefrontal cortex and a narrower 
region in the left prefrontal cortex, both of which were mainly located 
in the middle frontal gyrus. 



4. Discussion 

The present study investigated the prefrontal hemodynamic re- 
sponses to the forward and backward digit span tasks among veter- 
ans with PTSD and comorbid conditions, by using multichannel fNIRS. 
An event-related stimulation design was utilized to reveal the neural 
correlates of three working memory processes, namely the encoding, 
maintenance, and retrieval phases. By comparing with the results 
from the gender- and age-matched healthy controls, PTSD in veter- 
ans was found to be associated with distinct deactivations during the 
retrieval process. Significant prefrontal response differences between 
the two participant groups were seen predominantly in the right pre- 
frontal cortex as well as a narrow region in the left prefrontal cortex. 
The results of the present study and the possible neural mechanisms 
underlying these results are discussed in more detail below. 



4 A. Results from the healthy controls 



3.3. Correlations between behavioral and hemodynamic measures 

For both the forward and backward tasks, linear regression anal- 
yses between the individual behavioral and hemodynamic measures 
were conducted only in the retrieval phase. This was because while 
the behavioral measures were significantly different between the two 
participant groups, the hemodynamic measures showed significant 
group differences only in the retrieval phase, as described above. For 
the forward task, there was a positive, significant correlation (R = 0.42, 
p = 0.02) between the individual performance accuracies and the 
within-ROIs hemodynamic responses in the retrial phase, as shown 
in Fig. 9A. For the backward task, the correlation between these two 
measures was low and insignificant (R = 0.002, p = 0.99), which is 
shown in Fig. 9B. 



In the control group, the digit forward task evoked robust hemo- 
dynamic activations during the encoding and retrieval phases, which 
is consistent with a previous fMRI report using a similar stimulation 
protocol (Sun et al, 2005). The digit backward task evoked a much 
greater activation in the retrieval phase, which is also consistent with 
the previous fMRI report (Sun et al, 2005). Compared with the for- 
ward retrieval, the backward retrieval requires the digits to be held 
in short-term memory and a mental manipulation to be performed 
on them, and therefore is believed to engage more executive pro- 
cesses. A greater hemodynamic activation in the backward retrieval 
is in agreement with a greater neural demand induced by this task. 

In the present study, the digit backward task also evoked a subtle 
peak of Hb02 changes in the maintenance phase, which, however, was 
not seen in the fMRI study of Sun et al. (2005). This may be presumably 
because the present study had a longer duration of maintenance (10s) 
than in the previous study (6 s). As a result, some participants in the 
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Fig. 4. Topographic images of the task-evoked prefrontal activations (derived from A[Hb0 2 ]) in the control group. First row: mean prefrontal activations (/3-maps, in microM) 
evoked by the digit forward (DF) task in the encoding, maintenance, and retrieval phases. Second row: mean prefrontal activations (/3-maps, in microM) evoked by the digit 
backward (DB) task in the encoding, maintenance, and retrieval phases. Third row: mean differences of activations (/3-maps, in microM) between the backward and forward tasks 
(DB-DF) in the encoding, maintenance, and retrieval phases. Fourth row: t-statistical maps (t-maps) which outline the regions of significant difference (paired t-test, p < 0.05) 
between the backward and forward tasks in the encoding, maintenance, and retrieval phases. 



present study might have started to reverse the digit string in the 
maintenance phase, which could evoke higher Hb0 2 changes. 

The healthy controls also showed significant task differences in 
behavioral measures, as their accuracies in the backward task were 
significantly lower than those in the forward task. The reduced accu- 
racies in the backward task are in agreement with the more complex 
processes involved in this task, as described above. 



evolutions and magnitudes. The veterans with PTSD also showed sim- 
ilar performance accuracies across the two tasks. It seems that there 
was a ceiling effect in the veterans with PTSD when completing the 
forward and backward tasks: while their performance was already 
poor in the forward task compared with their healthy counterparts, it 
did not become worse in the backward task. This ceiling effect could 
also explain the high similarity of the hemodynamic responses to the 
two tasks in veterans with PTSD. 



4.2. Results from the veterans with PTSD 

For the veterans with PTSD, both the forward and backward tasks 
evoked hemodynamic activations during the encoding and mainte- 
nance phases, which were consistent with the healthy controls, fol- 
lowed by distinct deactivations in the retrieval phase. The hemody- 
namic responses to the two tasks were similar in terms of temporal 



4.3. Differences between the two participant groups 

In the present study, statistical comparisons of task-evoked hemo- 
dynamic responses between the two participant groups revealed that 
veterans with PTSD had significantly abnormal prefrontal responses 
compared with healthy controls while processing the working mem- 
ory tasks. A remarkable finding is that the healthy controls had nearly 
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Fig. 5. Mean task-evoked prefrontal hemodynamic responses (channel 30 in Fig. 1 A) in the PTSD group: (A) mean hemodynamic responses to the digit forward task, in which the 
mean Hb0 2 changes show a predominant peak in the maintenance phase followed by a distinct dip in the retrieval phase; (B) mean hemodynamic responses to the digit backward 
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baseline, and the gray bars in the bottom indicate the three phases of the task. It is noted that the actual retrieval time was slightly variable among participants. Therefore the 
retrieval phases labeled in both panels are schematic and approximate. 



equal activations in the two hemispheres during all three phases 
(see Fig. 3), but veterans with PTSD showed dominant activations 
in the left DLPFC in the encoding and maintenance phases (see Fig. 5). 
These observations were consistent with a study by Schiffer et al. 
(1995) that demonstrated differences in activity between right and 
left hemispheres evoked by negative vs. neutral (i.e., not trauma- 
related) memories in patients with a history of trauma, but not in 
controls. According to their hemispheric emotional valence (HEV) 
model, the right frontal cortex is more associated with the processing 
of negative emotions, and the left frontal cortex is more associated 
with the processing of positive and neutral emotions. When recalling 
emotional memories, this right-left difference was significant in pa- 
tients with a history of trauma, whereas control subjects had smaller, 
insignificant changes (Schiffer et al, 1995). The digit span tasks used 
in our study represented a neutral (i.e., not trauma-related) memory, 
and should therefore show more lateralized activations to the left pre- 
frontal cortex of PTSD participants than in controls, as predicted by 
the HEV model. Although another study (Schiffer et al, 2007) pointed 
out that HEV-laterality is more individually determined, one could 
still expect a majority of PTSD patients to have predominantly left- 
sided prefrontal activations during a neutral working memory task, 
as we found in our study. 

More importantly, veteran PTSD was found to be associated with 
distinct deactivations during the retrieval process, especially on the 
right prefrontal cortex (see Figs. 6 and 7). The identified ROIs that 
showed significant group differences included the right prefrontal 
cortex, localized at the right middle frontal gyrus, and a narrow re- 
gion in the left prefrontal cortex, localized within the superior frontal 
gyrus in the forward task and within the middle frontal gyrus in the 
backward task. While the present study could not determine the exact 
neural mechanisms underlying this deactivation involved in memory 
retrieval, several possible interpretations are discussed below. 

First, the veterans with PTSD had significantly poorer performance 
accuracies, indicating that they might have greater difficulty than 
the healthy controls in completing both the forward and backward 
tasks. However, the hemodynamic deactivations during retrieval in 
veterans with PTSD are unlikely a result of greater task difficulty. 
Previous studies (Hoshi et al, 2000; Gerton et al, 2004) on a group 
of healthy controls had demonstrated that a greater task difficulty 
(or with a longer length of digit string) would not cause deactivation 
during the retrieval process. To further confirm this conclusion, we 



also retested nine healthy participants in the control group with eight- 
digit forward tasks. When the task difficulty was increased from six 
digits to eight digits, these participants still showed clear activations 
(significant Hb0 2 increase) in the retrieval phase (Fig. 10). 

Second, as reviewed in the Introduction, a few neuroimaging stud- 
ies of working memory have consistently found that hypoactivation 
in the prefrontal cortex is associated with PTSD (Matsuo et al, 2003; 
Clark et al, 2003; Moores et al, 2008). These studies used block de- 
signs to measure the accumulated neural correlates of the working 
memory tasks. In the present study, if we assume that the accumu- 
lated neural correlates in a block design could be estimated as a sum- 
mation of event-related changes across all of the transit memory 
processes (encoding, maintenance, and retrieval), the hemodynamic 
activations during encoding and maintenance would be reduced by 
the deactivation during retrieval in veterans with PTSD, that should 
lead to results consistent with the previous neuroimaging studies. 
Thus, some of the findings in the previous studies can be applied to 
interpret the results in this study. The study by Matsuo et al. (2003) 
using the same fNIRS technology found a strong correlation between 
the prefrontal memory-related hypoactivation with the lower scores 
on attention and concentration among PTSD patients (victims of the 
Tokyo Subway Sarin Attack). They concluded that the prefrontal hy- 
poactivation might be secondary to a reduced attention and concen- 
tration capacity in PTSD patients; namely, the prefrontal cortex was 
not efficiently concentrated or engaged in the given memory task. The 
results in the present study revealed, however, that the prefrontal 
'hypoactivation' could be a composite outcome of activations during 
the encoding and maintenance phases followed by a deactivation dur- 
ing the retrieval phase. The deactivation cannot be explained simply 
as an insufficient engagement of the prefrontal cortex, because an 
insufficient brain activity or engagement theoretically should cause 
hemodynamic activation with reduced magnitude, rather than hemo- 
dynamic deactivation. 

Finally, we speculate that the distinct deactivation revealed by the 
present study may represent an active inhibition of prefrontal activi- 
ties in the retrieval phase (Tomasi et al, 2006). The PTSD participants, 
but not the controls, appeared to suppress prefrontal activity during 
memory retrieval. This deactivation was more pronounced in the right 
DLPFC. This finding is consistent with a previous fMRI study of PTSD 
patients performing a memory task of 1 2-word-pair associates (Geuze 
et al, 2008). Similar to our digit span working memory tasks, the 
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Fig. 6. Topographic images of the task-evoked prefrontal activations (derived from A[Hb0 2 ]) in the PTSD group. First row: mean prefrontal activations (/?-maps, in microM) evoked 
by the digit forward (DF) task in the encoding, maintenance, and retrieval phases. Second row: mean prefrontal activations (/3-maps, in microM) evoked by the digit backward 
(DB) task in the encoding, maintenance, and retrieval phases. Third row: mean differences of activations (/3-maps, in microM) between the backward and forward tasks (DB-DF) 
in the encoding, maintenance, and retrieval phases. Fourth row: t-statistical maps (t-maps) which outline the regions of significant difference (p < 0.05) between the backward 
and forward tasks in the encoding, maintenance, and retrieval phases. 



reported associative memory task consisted of neutral (not trauma 
related) material, and PTSD subjects also showed under-activation of 
the right prefrontal cortex during retrieval. Furthermore, both struc- 
tural and functional neuroimaging studies support that impairment of 
the right DLPFC is involved in PTSD, and that normalization of the right 
DLPFC is linked to recovery from trauma (Lyoo et al, 201 1 ). Repetitive 
transcranial magnetic stimulation (rTMS) of the left or right DLPFC is 
effective at reducing symptoms in PTSD (up to 3 months follow-up), 
but increasing cortical excitability in the right DLPFC induces a larger 
beneficial effect compared to the left (Boggio et al, 2010). Together 
with our findings, these studies suggest that normalization of DLPFC 
activity, particularly on the right hemisphere, might be a direction for 
developing treatments of PTSD. 



4.4. Correlations between the behavioral and hemodynamic measures 

For the forward task, a positive, significant correlation between 
the individual performance accuracies and hemodynamic responses 
in the retrieval phase was found for both groups. However, for the 
backward task, the correlation between these two measures was low. 
Such a low correlation in the backward task might be attributed to 
the potential ceiling effect in the PTSD group. As discussed above, 
although the backward task is apparently more difficult than the for- 
ward task, the veterans with PTSD showed small differences in both 
the behavioral and hemodynamic measures between the two tasks. 
In a scenario in which veterans with PTSD reach a response ceiling 
in the backward task, a low correlation between the behavioral and 
hemodynamic measures in this task is expected. A useful piece of in- 
formation obtained from the behavioral-hemodynamic correlations 
is that, for future neuroimaging studies, the digit forward task seems 
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Fig. 7. Topographic images (derived from A[Hb0 2 ]) showing the differences between the control and PTSD groups in response to the digit forward task. Top row: mean differences 
of activations (/3-maps, in microM) between the control and PTSD groups in the encoding, maintenance, and retrieval phases of the digit forward task. Bottom row: t-statistical 
maps (t-maps) which outline the regions of significant group differences (p < 0.01) in the encoding, maintenance, and retrieval phases of the digit forward task. 
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Fig. 8. Topographic images (derived from A[Hb0 2 ]) showing the differences between the control and PTSD groups in responses to the digit backward task. Top row: mean 
differences of activations (/?-maps, in microM) between the control and PTSD groups in the encoding, maintenance, and retrieval phases of the digit backward task. Bottom row: 
t-statistical maps (t-maps) which outline the regions of significant group differences (p < 0.01) in the encoding, maintenance, and retrieval phases of the digit backward task. 



a better task to use for objective assessment of the memory-related 
dysfunctions among PTSD patients. This task adapts well to partic- 
ipants with both normal and declined neurocognitive abilities. The 
task-evoked hemodynamic responses predict the participants' per- 
formance in a broad range. In contrast, the digit backward task seems 
too difficult for patients with impaired memory functions. 

4.5. Limitations of the present study 

We acknowledge several technical limitations in the present study. 
First, the fNIRS probe in the present study covered only a limited re- 
gion of the prefrontal cortex. The involvement of other cortical regions 
during the digit span tasks was not investigated. The present study 



can be improved by using a multi-patch probe to cover more cortical 
regions that deal with memory and attention. Second, in the GLM 
analysis, we used a standard HRF derived from BOLD fMRI (Glover, 
1999) as a surrogate to model the signals measured by fNIRS. This 
approach is reasonable because both fMRI and fNIRS detect similar 
hemodynamic changes in the brain. However, a previous study using 
a finger-tapping task did show that temporal profiles in Hb0 2 changes 
deviated from those of BOLD signals (Huppert et al, 2006). It means 
that an HRF derived from BOLD signals is not the best to model the 
Hb02 changes quantified by fNIRS. Lastly, the veterans with PTSD in 
the study had comorbid conditions, which may create hemodynamic 
signals that confound those resulting only from PTSD conditions. 
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Fig. 9. Correlations between the individual performance accuracies and hemodynamic responses to: (A) digit forward task, and (B) digit backward task. For each individual 
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Fig. 10. Mean prefrontal hemodynamic responses (channel 30 in Fig. 1A) to the eight- 
digit forward task among nine retested healthy controls. In this graph the solid lines 
represent the mean time courses of Hb0 2 , the dotted lines represent the mean time 
courses of Hb, the shaded regions represent the standard errors, the * symbols indicate 
the period of significant changes (p < 0.01) from the baseline, and the gray bars in 
the bottom approximately indicate the three phases of the task. The mean time course 
of HbC>2 shows two robust peaks in the encoding and retrieval phases, similar to 
the responses evoked by the six-digit forward task (Fig. 3 A). It demonstrates that for 
healthy controls, an increased task difficulty (from six digits to eight digits) still causes 
prefrontal activation rather than any deactivation. 



among the veterans with PTSD. The present study demonstrated that 
fNIRS could be a portable and complementary neuroimaging tool to 
study the cognitive dysfunctions associated with PTSD. We expect 
that a larger coverage of the cortex with more optodes and measure- 
ment channels will permit us to simultaneously reveal essential cor- 
tical regions involved in cognitive processes. In the meantime, future 
studies should seek a better HRF derived from fNIRS signals directly, 
while much effort on exclusion of subjects with comorbid conditions 
from the study is needed in order to accurately investigate prefrontal 
responses to memory phases in subjects with PTSD only. 
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5. Conclusions and future work 

In the present study, multichannel fNIRS was utilized to assess the 
involvement of the prefrontal cortex in working memory processes 
among veterans with PTSD and age-/gender-matched healthy con- 
trols. The veterans with PTSD, but not the healthy controls, were found 
to have distinct deactivations during the retrieval phase of working 
memory, which might implicate an active inhibition of prefrontal 
neural activity. The deactivation in the retrieval phase was more pro- 
nounced in the right dorsolateral prefrontal cortex than the left side 
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